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A class of fc-particle observables in a two-mode system of Bose particles is characterized by typical¬ 
ity: if the state of the system is sampled out of a suitable ensemble, an experimental measurement of 
that observable yields (almost) always the same result. We investigate the general features of typical 
observables, the criteria to determine typicality and finally focus on the case of density correlation 
functions, which are related to spatial distribution of particles and interference. 
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I. INTRODUCTION 

The experimental realization of Bose-Einstein conden¬ 
sates raised great theoretical interest. A system of Bose 
particles in one or few single-particle states (modes) is an 
important workbench for fundamental concepts in quan¬ 
tum mechanics and statistical physics [THU. A large num¬ 
ber of particles distributed among two different modes, 
for example, enables one to perform a full quantum 
double-slit experiment in a single experimental run [5]. 
Moreover, it has been proposed that fluctuations in the 
interference patterns can probe interesting characteristics 
of many-body systems ims]. 

Interference is an interesting example of a property 
that weakly depends on the choice of the state of the sys¬ 
tem. In a two-mode system, second-order-interference 
properties are similar as far as one considers a num¬ 
ber state or a phase state 0 H HO], while first-order 
properties are very different. These features explain 
why an interference pattern can be experimentally ob¬ 
served in single experimental runs by measuring the par¬ 
ticles’ positions [HHnunMa, although in the case of 
number states the offset of the pattern fluctuates ran¬ 
domly, so that averaging over a few experimental runs 
yields a flat density profile. It Is also Interesting to note 
that, when the number of particles Is large, there are 
interference-related observables whose experimental mea¬ 
surement yields the same result at each run with over¬ 
whelming probability PlHIHl- 

These considerations lead us to a general definition of 
typicality of an observable. Typicality is a mathematical 
concept related to the phenomenon of measure concen¬ 
tration m- It has been used in many emerging phe¬ 
nomena in physics and other sciences, with interesting 
applications on the structure of entanglement In large 
quantum systems |28H32j . and the search for a quantum 
mechanical justification of some primary statistical me¬ 
chanical concepts [HHl - HH] . In this article, we shall apply 
the notion of typicality to a two-mode Bose system with 
a fixed number of particles. We will study the proper¬ 
ties of observables with respect to uniform sampling of 
a suitable Hilbert subspace. An observable will be de¬ 
fined typical if each experimental run, performed on any 
state of the subspace, would yield the same result with 


overwhelming probability. According to the probabilis¬ 
tic interpretation of quantum mechanics, the expectation 
value of an observable provides information on the aver¬ 
age of experimental runs on the (pure or mixed) state. 
We shall build on the results of Ref. [55] and prove that 
If the observable is typical, there exists an expectation 
value that contains information on (almost) each single 
experimental run, rather than on the average result. 

This article Is organized as follows. In Section |TT| we 
specify the properties of the statistical ensemble and In¬ 
troduce the general formalism, which leads to the defi¬ 
nition of typicality of an observable. In Section |III| an 
analysis of the structure of fluctuations for a general k- 
particle observable is performed. Quantitative criteria 
will be given to determine typicality and obtain Informa¬ 
tions on higher-order fluctuations. Section |IV| includes 
an application of typicality criteria to density correla¬ 
tion operators, with the case study of counterpropagating 
plane-wave modes and an extension to expanding modes 
in the far-held regime. 

II. TYPICALITY OF AN OBSERVABLE 

Let us consider a system of bosons, without internal 
degrees of freedom. The system can be described in a 
second-quantization picture by introducing the annihila¬ 
tion and creation Held operators ^(r) and ^^(r), satis¬ 
fying canonical equal-time commutation relations: 

['I'(r),'k(rO] = 0, (1) 

[^(r),^^(r')] = 5{r-r'). (2) 

Once an orthonormal basis for the single-particle Hilbert 
space has been chosen, the Held operators can be ex¬ 
panded as sums of mode operators, which create or an¬ 
nihilate a particle in one of the basis states. We are 
going to analyze a system made up of N bosons, dis¬ 
tributed among two orthogonal modes, a and b, with 
single-particle wave functions and ipbif)- A useful 

basis for the description of the Y-particle Hilbert space 
Is formed by the Fock states 
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in which the two modes have well-dehned occupation 
numbers. We are assuming that N is even for simplic¬ 
ity, but this specihcation is immaterial in the large-A^ 
limit. In a second-quantized formalism, Fock states are 
obtained by applying a sequence of mode creation oper¬ 
ators to the vacuum |f2): 


K) 


1 

^{N/2 + i)\{N/2-iy. 


(~t)Af/2-h^(5t)W/2-^|^)^ 


( 4 ) 

Since ipaif) and ipbi'i') are orthonormal, the mode oper¬ 
ators 


a = 


dril}l{r)^{r), 


h = 


dr (5) 


satisfy the canonical commutation relations 

[d,ay] = [b,P] = l, (6) 

and all the operators of mode a commute with those of 
mode b. The number operators Na = 0 ^ 0 , and N}, = w b 
count the numbers of particles in each mode. 

We are interested in the properties of an ensemble of 
states randomly sampled from the n-dimensional sub¬ 
space 


■Hn = span{|£), |£| < n/2}, (7) 

spanned by the Fock states with 0 < n < N + 1 {n is odd 
for simplicity.) The case n = 1 represents a degenerate 
ensemble in which only the state |£ = 0) has nonvanishing 
probability [UlMMlIll]. We shall consider the large-n 
case [5^. If a BEC made up on N particles is “evenly” 
split between the two modes, one expects n = 0{'\fN). 
However, more general cases are possible [26]. In the 
following we will assume that 

n = o{N). (8) 

The assumption of uniform sampling is clearly a simpli¬ 
fying one: the number of states that are involved in the 
description and their amplitudes depends on the experi¬ 
mental procedure, yielding the separation of the conden¬ 
sate in the two modes [IS]. However, it will emerge that 
our main results are qualitatively unchanged for a large 
and relevant class of probability distributions on Hn- 
The average of the projection on a random vector state 
|$ 7 v) G din uniformly sampled on Hn yields the (micro- 
canonical) density matrix 

Pn = |‘f’Ar)($Ar| = - ^ |^)(£| =: -P„, (9) 

n ^' n 

h|<n/2 

where Pn is the projection onto the subspace Hn- In the 
following, only the density matrix (|^ will explicitly ap¬ 
pear in our analysis and calculation. The definition and 
use of the random state \^n) in Eq. is superfluous and 
can be dispensed with, in accord with the prescription of 


Ockham’s razor. In this respect, it is worth stressing 
that no hypothesis of decoherence will be made and if 
one wants one can safely assume that in each experimen¬ 
tal run a wave function describes the condensate, which 
is therefore in a pure state. 

Given an observable A, the statistical average of its 
expectation value reads 


A:=Tr{pnA) = - V {e\A\i). 
n ^' 

\e\<n /2 


Its quantum variance is 


:= Tv{pnA^) 


TT{pnA) 
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( 10 ) 


( 11 ) 


The signihcance of this quantity in the context of ensem¬ 
ble statistics becomes clear once it is decomposed in the 
sum of two contributions [251 SS] : it contains both the 
classical and quantum uncertainties of the observable A 
in the microcanonical state Therefore, the asymp¬ 
totic condition 


5A = o{A), (12) 

for N ^ oo, ensures that in the overwhelming majority of 
cases the experimental measurement of the observable A 
will fluctuate within an extremely narrow range around 
the average expectation value A. Thus, the outcome of 
a measurement of A is almost always the same (and it 
equals its average) for every experimental run in the en¬ 
semble. We call this property typicality of the observable. 

In the following sections we will define and characterize 
conditions for a fc-particle observable to be typical, and 
analyze in detail the case of spatial correlation functions, 
which are related to interference. 


III. CONTROL OF FLUCTUATIONS FOR A 
fe-PARTICLE OPERATOR 


In a second quantization formalism, the held operators 
^ and ^1 can be used to build up many-body observables 
m- The simplest ones are the Hermitian 2fc-point func¬ 
tions 


Gkiri, ..., rfc) := 4'l(r,)^(ri). 


( 13 ) 


Following Eq. (10), one can observe that the ensemble 


average of Gk can be computed for any n once its expec¬ 
tation value over the A-particle-two-mode Fock states \t) 
is known. Thus, when the held operators in (131 are ex¬ 


panded in orthogonal modes, only the terms formed by 
operators associated to modes a and b are relevant in Gk- 
Moreover, operator products with a different number of 
al’s and d’s (respectively d^’s and 5’s) give vanishing con¬ 
tributions to {£\Gk\'i)- The relevant terms in (13) can be 
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reordered to be expressed as number operators Na,b- The 
average eventually reads [we will assume k = 0(1)] 

k 

Gkiri,...,rk) = Fm{ri,... ,rk) 


m—0 


71 ^ ^ \ 2 

\i\<Ti/2 A=0 


k — m— 1 




B=0 


with 


Fm{ri, ... ,rfe) = |$m(ri,... ,rfe) 


(14) 


(15) 


where $m(ri,..., r^) is, apart from a normalization fac¬ 
tor, the symmetrized fc-body wave function with m par¬ 
ticles in mode a and k — m particles in mode b: 

^m{ri, ...,rk) = 

^ ^ '0a(^CT(l)) ■ ■ • '*Aa(^cr(m))'*Ab(^fT(m+l)) • ■ ■ '4^b{f'a-{k)) i 

(7 

(16) 

a denoting permutation of k elements. 

An important class of A:-particle observables can be 
obtained by integrating the 2_)-point functions, for j < fc, 
with a multiplicative kernel Aj m- 

k p 

Ak = Yj dri...drjAj{ri,...rj)Gj{ri,...rj). (17) 
3=0 d 

The ensemble average of A can be immediately com¬ 
puted inserting the general results (14). Thus, since 
Gj = 0{N^) when k = 0(1), we have 


Ak= dri...drkAkiri,...rk)Gkiri,...rj)+0{N^ ^). 

(18) 

In order to study the behavior of the variance 0 and 
determine whether Ak is typical or not, we should com¬ 
pute its square, which involves at the highest order in N 
a product of Gk functions, that can be recast into a single 
4fc-point function by normal-ordering the field operators: 

Al = J dri...dr2k{Ak{ri,...,rk)Ak{rk+i,---,r2k) 

xG2fe(ri,...,r2fe))+0(iV2fe-i). (ig) 

From 0 -([T^ we get an expression for the variance 

SAl = J dri .. .dr2k{Akiri,... ,rk)Akirk+i, ■ ■ ■,r2k) 

X7fc(r'i,...,r-2fe))+0(iV2'=-i), (20) 

with 


lk{ri,---,r2k) ■= G2fe(ri,...,r2fe) 


- Gfc(ri,... ,rfc) Gfe(rfc+i,...,r2fc), 

( 21 ) 


that depends only on the ensemble averages of Gk and 
G 2 k- Notice that the order comes both from con¬ 

tributions in ( |17[ ) with j < k and from normal ordering. 

Our aim is to find whether, for some choice of the sam¬ 
pled Hilbert subspace Hn , the standard deviation (201 is 


of smaller order with respect to the average Ak in the 
asymptotic N ^ oo regime. This behavior is in ac¬ 
cord with the definition of typicality of the observable 
Ak given in Sec.|^ First, let us observe that, since Gk is 


polynomial in N and n, inserting the general result (14) 
into SAl yields a polynomial function of degree 2k in the 
number of particles and the dimension of the sampled 
subspace: 


fc 2k—p 


= E E A7 


p—0 q—0 




( 22 ) 


Due to the symmetry of the summand around i = 0, 
only terms with even q are present in (22). It 


is evident that, as far as n = o{N), it is necessary and 
sufficient for Ak to be typical that 


D 


Ak) 
2k,0 


j dri .. .dr2kAk{ri,... ,rk)Ak{rk+i, ■ ■ ■,r2k) 

k 

Y FM{ri,...,r2k)- 

^ ^ ^m' (’^li ■ • ■ • ■ • ■) 'f^2k) j (2^) 

'm '—n X 




vanish, namely. 


J^Ak) _ Q 

^2k,ei ~ '0- 


(24) 


In this case, since Ak = 0{N^), the relative fluctuations 
are 


SAk 

Ak 




(25) 


Fluctuations that scale like related to normal or¬ 

dering and to the very definition of Ak, are ensemble- 
independent, in the sense that they are present even in 
degenerate distributions of states. Linear fluctuations in 
n are clearly related to the dimension of the sampled sub¬ 
space, and therefore strongly depend on the definition of 
the ensemble. It is interesting to note that, if n = 0{N°‘), 
two qualitative regimes can be distinguished: i) when 
a < 1/2, the relative fluctuations scale like ii) 

while if a > 1/2, they asymptotically vanish like 
i.e. more slowly. This reasoning is based on the assump¬ 
tion that D 2 Y 2 2 does not vanish: we will see in the 
following how this condition can be checked and discuss 
in the next section a relevant class of exceptions. 

Due to the form of the symmetrized products of wave 
functions (16), the typicality condition ([24| can be recast 


in a more convenient form: in particular, we can dispose 
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of the integral over 2k variables by observing that the 
functions , J' 2 fe) can be decomposed as 

$M(r'i,.. ■,r2k) = 
mm{M,k} 

■ ■ • T '^2k'} • 

m—nmx{0,M—k} 

(26) 

This result is related to the cluster decomposition princi¬ 
ple in Bose systems [H] . When one computes the square 
modulus Fm = |4>Mp of (26), that enters the variance in 


(21), the term ^M-mFm appears. This contribution 


exactly cancels with the subtracted terms in (21), since 


it can be obtained by a change of summation indices 

k 

Y .. . ,rk)Fm{rk+i, ■. ■,r2k) = 

k M 

EE FM-m{ri, ..., rk)Fm{rk+i, • • •, r2k) 

M—0 m—0 
2k k 

+ E E FM-m{ri ,..., rk)F^{rk+i ,..., rafc). 

M—k-\-l m—M—k 

(27) 

Accordingly, a compact form of the typicality condi¬ 
tion (24) reads 


r{Ak) 


n-iAk) _ n 

rri' 


2k min{M,k} 

jjiAk) _ j-y 

'^2fc,0 / ^ ^ "^M—m,M— 

M—0 m'7^m=max{0,M—fe} 

(28) 

where the coefficients Irn,m' are integrals over k position 
variables of products of the type 


A 


= fdri... drkAkiri, rfc)($;,,4>„) (n,..., rj,). 

(29) 

Given the integral kernel Ak-, which determines the high¬ 
est order in N of the observable Ak, and the mode wave 
functions ipa,bi'i')y h is sufficient to compute the integrals 
(29) to check whether the observable is typical. 


It is also interesting to analyze the structure of 
D 2 Y 2 2 > view of the following discussion on spa¬ 

tial correlation function. If 7 I 2 E 2 2 nonvanishing, 
the transition between an ensemble-independent 
and an ensemble-dependent behavior of fluctuations oc¬ 
curs for n = 0 {N^A'j^ jn order to analyze this quan¬ 
tity, one should take into account terms of order and 
Nk- 2^2 general expression (14) 


k 

Gfc(ri,...,rfc) = Y Fm{ri,...,rk) 
m—0 

Y) ' \ 2 J 24 


X 


Ar\ k—2 2 

N\ 


(k^ — k{4m -I- 1) -I- 4m^) 


(30) 


to be used in the computation of (21). Integration over 


the kernel Ak and application of the same change of in¬ 
dices leading to Eq. ^7\ yield the result 


D 




{Ak) 

2fc—2,2 ^2 

2k min{M,fc} 


+ E E i^M- 

M—O m'^m—max{0,M—fc} 

X (2fc2 - fc(4M -f 1) -f 2M2)] |, 


m,M—m' m,m' 


(31) 


where the X integrals have been defined in (29), and 


= J dri... drkAkiri, ...,rk) 


X E “ 2m). (32) 


m—0 


The summation appearing in (32) contains an equal num¬ 


ber of positive and negative term. Thus, it can vanish 
if the mode structure is properly chosen. In particular, 
it vanishes when the E^’s are invariant with respect to 
exchange of the two mode wave functions: 


Frr, = Fl. 


k—m 


Y ■^™(^ “ 2m) = 0. 


(33) 


m—0 


This condition is always valid in a “double-slit” BEG in¬ 
terference experiment [SHMSS], where the two modes 
are identically prepared (within experimental accuracy) 
and then let to interfere. It is easy to check that the 
above condition is satisfied whenever |V^a(j')| = \'4’bi'f')\ 
at all points. In turn, this is a consequence of the invari¬ 
ance of Fm under local phase transformations ipa,bii') -A 
F^F)^a,bir). 


IV. TYPICALITY OF DENSITY 
CORRELATIONS 

In this section we will specialize the general results 
obtained for an observable of the form 0 to a partic¬ 
ular class, related to the spatial interference of conden¬ 
sates. This is usually accessible to experimentalists and 
provides an interesting example of typical behavior. We 
shall analyze the integrated density correlation functions 

Ckixi,... ,Xk-i) = j drpir)pir + xi)... pir + Xk-i) 
= J drGkir,r + Xi,...,r + Xk-i) 

-tO(A'=-i), (34) 


where /5(r) = Gi(r). In agreement with Eq. (17), it 


is clear from the normal-ordered form that the highest- 
order integral kernel for this class of observables is 


k-l 


Cfe(ri, ...,rk)= (5(n+i - Ti - Xi). (35) 
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The singularities arising in the normal ordering of (34) of the modes provides interesting information als o o n the 


can be avoided by smearing all densities around the 
points (r + Xi) with functions that take into account 
the finite experimental spatial resolution. We will focus 
on two modes that are paradigmatic in the description 
of interference of Bose-Eisntein condensates, namely two 
counterpropagating plane waves 


V'a(r') = 


(36) 


and eventually extend the results to the case of two 
modes that are spatially separated at the initial time, 
and are let to expand and overlap: this is a more realis¬ 
tic description of the experiments. 

In order to satisfy the typicality condition (281, it is 


generally not necessary that the X integrals themselves 
vanish. However, in this case we can verify the stronger 
condition that all the integrals with m! ^ m are iden¬ 
tically zero. Let us compute the general form of these 
integrals by using (351 


^^3% = J dn--.drfeCfc(ri,...,rfc)($);,$„)(ri,...,rfc) 

= J r + xi,...,r + Xk-i), (37) 

and consider for definiteness m' > m. In this case, the 
function ..., r^) is the sum of products of k 

mode wave functions and k complex conjugates, with the 
number of V’a’s exceeding the number of tpaS hym' — m. 
The structure of the products reads 

(T—l T—1 

Z X 

X n (rjc ) n ) (38) 

C=i ?=i 

with S + T + Z + X = k, S + Z = m' and S + X = m, 
which implies 


Z — X = m! — m 


(39) 


for all products in For the plane-wave modes 

(36), IV'a.bl = 1 and Inserting these 

results in the general form of the products (38) and inte¬ 
grating over r as in (37) yields 


z X . 

JJ" g2ifco / ^^g-2i(m'-m)fco r _ q 

C=1 «=1 

(40) 


Thus, all the contributions to the integral (37) identically 


vanish for every k and every set of points ixi, ..., X}~_i). 
This implies that condition 


n(Cfc) _ 

^2fc.O “ ^ 


(41) 


is satisfied by all density correlation functions of the form 
(34), which are thus typical for n = o{N). The structure 


N'i.k 2^2 fluctuations, arising from (|21[) and 


the general expression (30). Since |'0a| and |'(/'b| are iden¬ 
tically equal to one, the mode wave functions satisfy the 
symmetry condition (33) for the functions. This im¬ 


plies that the integrals, defined as in Eq. (32), in- 


dentically vanish, which, together with the cancellation 
of the I integrals, leads to the result 


Jj{Ck) _ Q 
^2k-2,2 — 


(42) 


Thus, the highest order of ensemble-dependent contri¬ 
butions to the variance SC'^ is indeed which 

should be compared with the order of the 

ensemble-independent fluctuations. This means that if 
n = 0{N°‘), the relative fluctuations behave like 


5Cl = 


O(^) for a <3/4, 

O = C>(iV4(“-i)) for a > 3/4. 


(43) 

The transition between an ensemble-independent regime 
of relative fluctuations to an ensemble-dependent one 
does not take place, as one would expect, at n ^ Vx, 
but extends up to n ~ This behavior, which is 

general for density correlation functions, generalizes the 
results obtained in |26) for the Fourier transform of the 
second order density correlation function 


C 2 {x) = J drp(r)p{r + x). 


(44) 


We include a brief discussion on this operator to clarify 
what implications typicality has from an experimental 
point of view. The expectation value of (44) is generally 
given by 

C 2 {x) = ^ X! / drF^{r,r + x) + 0{N) + Oin"^), 

(45) 

which, in the case of plane waves, specializes to 


1 


C 2 {x) ~ iV^ 1-1-4 cos(2fco • x) 


(46) 


Due to typicality, the function (46) represents the over¬ 


whelmingly probable experimental result of a measure¬ 
ment of the observable < 72 ( 3 ;), unless n = 0{N). This 
function is the two-point density correlation of a classi¬ 
cal density (see discussion in Ref. [26) 1 


p{x) = 2N‘^ cos^(fco ■ X + (j)). 


(47) 


Among all the parameters that determine the typical 
experimental outcome of a density measurement, one, 
namely the offset (/ of the interference pattern, is not 
determined by typicality, since correlation functions do 
not depend on it. In Figure a possible outcome of a 
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p{x) 



FIG. 1: Outcome of a density measurement in a two plane- 
wave mode system, in the asymptotic regime. The visibility 
and the period of the interference pattern are (almost) fixed 
by typicality, while the offset of the pattern, and thus the 
value of density at the origin, fluctuates randomly. The order 
y/N of the visibility fluctuations refers to the case n = 0{N°‘) 
with a < 3/4. 


density measurement in a two plane-wave mode system 
is represented, with its typical features highlighted. 

Let us finally discuss the typicality properties of den¬ 
sity correlation functions for a two-mode system that is 
closer to actual experimental implementations, and find 


out in which cases the results (41)-(421 can be general¬ 
ized. Let us consider two spatially separated modes pa(z) 
and pb{z) in one dimension, which are concentrated re¬ 
spectively around positions Za and z;,, their initial width 
being much smaller than their distance. The spatial sep¬ 
aration implies that the convolution 


dz''lp*aiz')pb{z' ■ 


, 2 ) =: g.Fo.b(z)+iVabiz) 


(48) 


with Fab and (pab real functions of z, is peaked around 
a point z = zq, while the convolutions of p* with pa 
and of V'b with pb are peaked around z = 0. If the 
trapping potential is turned off at the initial time t = 0, 
the particles evolve under the free Hamiltonian F[q = 
—(P/dz'^. In absence of collisions, the A:-point functions 
at t > 0 can be obtained by replacing 

Pa,b{^) pa,b{z,t) = ex.ip{-itHo)pa,biz). (49) 

We are interested in the large-time {far-field) regime, in 
which the evolved wave function are approximated by the 
asymptotic form 


lution (48) as 


P*a{k)PPk) = j (51) 

Using a quadratic approximation of Fab around its maxi¬ 
mum zo, the integral becomes Gaussian and the product 
PaPb reads 


P*a{zp)Pb{z,t) ^ c{zo,t)exp(-i — '\ 


2t / 
2 


X exp — 




(52) 


If the spatial period of the complex exponential is much 
shorter than the standard deviation of the Gaussian part, 
namely 


zo > 


27r 


IKbi^oW 


(53) 


the product (52) is consistent with the result for counter- 


propagating plane wave modes with time-dependent wave 
number fco(t) = zo/4t, modulated by a slowly-varying 
Gaussian envelope, whose width is of the same order 
as the standard deviation of the far-field-approxim ated 
Ipal"^ and IV'bP- Thus, the integrals of products (38) 


with the kernel (35) are exponentially suppressed like 
g-(zolU6l) . This means that when the relation (53) on 


the initial wave packets is satisfied, typicality condition 
(41) on correlation functions is fulfilled within a very 


good level of approximation. 

The generalization of the result (42), yielding the can¬ 
cellation of N‘^^~‘^rP terms in the variance, cannot be 
extended without further assumptions to the case of ex¬ 
panding modes in the far-field regime. In fact, the re¬ 
lation \pa{z,t)\ = \pb{z,t)\ might not even be approxi¬ 
mately verified in the large-time limit. Gondition 
which, together with the results discussed above, leads to 
^ 2 k -2 2 “ verified only if the low-momentum 

Fourier components of the two modes are approximately 
equal [see Eq. (50).] In the case analyzed in Ref. [26] , in 
which the correlation C 2 {z) has been studied for a sys¬ 
tem with two translated Gaussian modes, it turned out 

/ p-1 \ 

that L? 2 /c 2 — Oj expected from the present general 
discussion. 


V. CONCLUSIONS AND OUTLOOK 


pa,b{z,t) ~ 


A-Klt 


■^a,b(2j 


(50) 


where pa^b{k) = J dze~^’‘^pa,b(z) are the Fourier trans¬ 
forms of the initial modes. In the far-held approximation, 
the quantities PaPb their complex conjugates, which 
enter the typicality condition ( |^ through products of 
the type (38), are thus related to products of Fourier 


transforms, which can be expressed through the convo- 


In this article we have dehned and characterized the 
typicality of a generic observable in a two-mode Bose sys¬ 
tem. The results enable one to determine if an observable 
is typical once the mode wave functions are known, and 
to identify different regimes in the huctuations around 
the typical expectation value, as the dimension of the 
sampled subspace varies. The observable analyzed in the 
last section are experimentally accessible, and can pro¬ 
vide a test for typicality. 
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The identification of typicality criteria for observables 
helps one understanding which properties are shared by 
the vast majority of states and which ones have instead 
wide fluctuations. Remarkably, this distinction is central 
in determining “good” (macroscopic) observables in both 
classical and quantum statistical mechanics [55]. The 
relation between the results obtained in this article and 
statistical mechanics will be the object of future research. 

It would also be interesting to extend the formalism 
in order to include more general cases, like nonuniform 
samplings or randomly fluctuating modes. Other possi¬ 
ble avenues for future investigation will be the analysis of 
the dynamical effects of typicality [3^ and the statistical 
interpretation of recent experiments on phase randomiza¬ 


tion in condensates [3SH53], as well as the characteriza¬ 
tion of the typicality of entanglement in a Bose-Einstein 
condensate jSSj. 


Acknowledgments 

This work was partially supported by the Italian Na¬ 
tional Group of Mathematical Physics (GNFM-INdAM), 
and by PRIN 2010LLKJBX on “Gollective quantum phe¬ 
nomena: from strongly correlated systems to quantum 
simulators.” 


[1] A. S. Parkins and D. F. Walls, Phys. Rep. 303, 1 (1998). 

[2] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999). 

[3] A. J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 

[4] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

[5] L. Pitaevskii and S. Stringari, Bose-Einstein Condensa¬ 
tion (Oxford University Press, Oxford, 2003). 

[6] C. J. Pethick and H. Smith, Bose-Einstein Condensation 
m Dilute Gases, 2nd ed. (Cambridge University Press, 
Cambridge, 2008). 

[7] A. J. Leggett, Quantum Liquids: Bose Condensation and 
Cooper Pairing in Condensed-Matter Systems (Oxford 
University Press, Oxford, 2006). 

[8] M. R. Andrews, C. G. Townsend, H.-J. Miesner, D. S. 
Durfee, D. M. Kurn, and W. Ketterle, Science 275, 637 
(1997). 

[9] E. Altman, E. Dernier, and M. D. Lukin, Phys. Rev. A 
70, 013603 (2004). 

[10] A. Polkovnikov, E. Altman, and E. Dernier, Proc. Natl. 
Acad. Sci. USA 103, 6125 (2006). 

[11] V. Gritsev, E. Altman, E. Dernier, and A. Polkovnikov, 
Nature Phys. 2, 705 (2006). 

[12] A. Polkovnikov, Europhys. Lett. 78, 10006 (2007). 

[13] A. Imambekov, V. Gritsev, and E. Dernier, in Ultra-Cold 
Eermi Gases, Vol. 164 of International School of Physies 
“Enrico Fermi”, edited by M. Inguscio, W. Ketterle, and 
C. Salomon (lOS, Amsterdam, 2007), pp. 535-606. 

[14] V. Gritsev, E. Dernier, and A. Polkovnikov, Phys. Rev. 
A 78, 063624 (2008). 

[15] Z. Hadzibabic, S. Stock, B. Battelier, V. Bretin, and J. 
Dalibard, Phys. Rev. Lett. 93, 180403 (2004). 

[16] S. Stock, Z. Hadzibabic, B. Battelier, M. Cheneau, and 
J. Dalibard, Phys. Rev. Lett. 95, 190403 (2005). 

[17] Z. Hadzibabic, P. Kruger, M. Cheneau, B. Battelier, and 
J. Dalibard, Nature (London) 441, 1118 (2006). 

[18] P. Kruger, Z. Hadzibabic, and J. Dalibard, Phys. Rev. 
Lett. 99, 040402 (2007). 

[19] K. Goral, M. Gajda, and K. Rz§,zewski, Phys. Rev. A 
66, 051602 (R) (2002). 

[20] Y. Castin and J. Dalibard, Phys. Rev. A 55, 4330 (1997). 

[21] J. Javanainen and S. M. Yoo, Phys. Rev. Lett. 76, 161 
(1996). 

[22] J. 1. Cirac, C. W. Gardiner, M. Naraschewski, and P. 


Zoller, Phys. Rev. A 54, 3714(R) (1996). 

[23] T. Wong, M. J. Collett, and D. F. Walls, Phys. Rev. A 
54, 3718(R) (1996). 

[24] M. lazzi and K. Yuasa, Phys. Rev. A 83, 033611 (2011). 

[25] S. Ando, K. Yuasa, and M. lazzi, Int. J. Quant. Inf. 9, 
431 (2011). 

[26] P. Facchi, H. Nakazato, S. Pascazio, F. V. Pepe, and K. 
Yuasa, Phys. Rev. A 89, 063625 (2014). 

[27] M. Ledoux, The Concentration of Measure Phenomenon, 
vol. 89 of Mathematical Surveys and Monographs (Amer¬ 
ican Mathematical Society, 2001). 

[28] P. Hayden, D.W. Leung, and A. Winter, Commun. Math. 
Phys. 265, 95 (2006). 

[29] P. Facchi, U. Marzolino, G. Parisi, S. Pascazio and A. 
Scardicchio Phys. Rev. Lett. 101, 050502 (2008). 

[30] A. De Pasquale, P. Facchi, V. Giovannetti, G. Parisi, S. 
Pascazio, and A. Scardicchio, J. Phys. A: Math. Theor. 
45, 015308 (2012). 

[31] P. Facchi, G. Florio, G. Parisi, S. Pascazio, and K. Yuasa, 
Phys. Rev. A 87, 052324 (2013). 

[32] F. D. Cunden, P. Facchi, and G. Florio, J. Phys. A: Math. 
Theor. 46, 315306 (2013). 

[33] S. Lloyd, Black Holes, Demons and the Loss of Co¬ 
herence: How Complex Systems Get Information, and 
What They Do with It, PhD thesis. Rockefeller Univer¬ 
sity (1991). 

[34] H. Tasaki, Phys. Rev. Lett. 80, 1373 (1998). 

[35] J. Gemmer and G. Mahler, Eur. Phys. J. B 31, 249 
(2003). 

[36] S. Popescu, A. J. Short, and A. Winter, The foundations 
of Statistical Mechanics from Entanglement: Individual 
States vs. Averages, quant-ph/0511225 (2005). 

[37] S. Popescu, A. J. Short, and A. Winter, Nature Phys. 2, 
754 (2006). 

[38] S. Goldstein, J. L. Lebowitz, R. Tumulka, and N. Zanghi, 
Phys. Rev. Lett. 96, 050403 (2006). 

[39] A. Sugita, Nonlinear Phenom. Complex Syst. 10, 192 
(2007) cond-mat/0602625 (2006)]. 

]40] P. Reimann, Phys. Rev. Lett. 99, 160404 (2007); 101, 
190403 (2008). 

[41] M. Rigol, V. Dunjko, and M. Olshanii, Nature (London) 
452, 854 (2008). 

[42] J. Cho and M. S. Kim, arXiv:0911.2110 (2009). 

[43] S. Sugiura and A. Shimizu, Phys. Rev. Lett. 108, 240401 




(2012); 111, 010401 (2013). 

[44] G. S. Paraoanu, J. Low. Temp. Phys. 153, 285 (2008); 
Phys. Rev. A 77, 041605 (R) (2008). 

[45] K. Maussang, G. E. Marti, T. Schneider, P. Treutlein, Y. 
Li, A. Sinatra, R. Long, J. Esteve, and J. Reichel, Phys. 
Rev. Lett. 105, 080403 (2010). 

[46] P. Facchi, H. Nakazato, S. Pascazio, F. V. Pepe, G. A. 
Sekh, and K. Yuasa “Phase randomization and typicality 
in the interference of two condensates” arXiv: 1410.1381 
[quant-ph]. 

[47] A.L. Fetter, and J.D. Walecka, Quantum Theory of 
Many-particle Systems (McGraw-Hill, New York, 1971). 

[48] S. Weinberg, The Quantum Theory of Fields, Volume 1 
(Gambridge University Press, Gambridge, 1996). 

[49] T. Schumm, S. Hofferberth, L. M. Andersson, S. Wilder- 
muth, S. Groth, 1. Bar-Joseph, J. Schmiedmayer, and P. 
Kruger, Nature Phys. 1, 57 (2005). 

[50] S. Hofferberth, 1. Lesanovsky, B. Fischer, J. Verdu, and 
J. Schmiedmayer, Nature Phys. 2, 710 (2006). 

[51] S. Hofferberth, 1. Lesanovsky, B. Fischer, T. Schumm, 


and J. Schmiedmayer, Nature (London) 449, 324 (2007). 

[52] S. Hofferberth, 1. Lesanovsky, T. Schumm, A. Imam- 
bekov, V. Gritsev, E. Dernier, and J. Schmiedmayer, Na¬ 
ture Phys. 4, 489 (2008). 

[53] T. Langen, R. Geiger, M. Kuhnert, B. Rauer, and J. 
Schmiedmayer, Nature Phys. 9, 640 (2013). 

[54] G. Teschl, Mathematical Methods in Quantum Meehanics 
(American Mathematical Society, Providence, 2009). 

[55] L. D. Landau, and E. M. Lifshitz, Statistical Physics, 
Part 1 (Pergamon Press, Oxford, 1980). 

[56] F. Benatti, R. Floreanini, and U. Marzolino, Phys. Rev. 
A 85, 042329 (2012). 

[57] We shall focus on this case for simplicity, but the formal¬ 
ism, as well as the results, can be generalized to kernels 
involving derivatives. On the other hand, momentum- 
dependent observables can be more easily treated by 
using correlation functions in momentum space, con¬ 
structed with the Fourier-transformed field operators 
l'(p) = /dre-‘P”4'(r). 


